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CHAPTER I 


INTRODUCTION 


In this study, an analysis of electromagnetic backscatteri ng from 
an open ended inlet cavity is performed, and the calculations based on 
this analysis are compared with experimental measurements. 

The wave scattered by the inlet cavity when it is illuminated by an 
external source is composed mainly of two effects. The first is due to 
the diffraction from the rim at the open-end of the inlet and the second 
is due to the coupling of the incident electromagnetic energy into 
cavity modes which bounce around inside the cavity before radiating from 
the open end. The inlet cavity considered here is assumed to have 
perfectly conducting walls. 

The analysis is done on an inlet model with overall dimensions of 
several wavelengths; therefore, there are quite a few propagating modes 
inside the cavity. The modal analysis is complemented with the high 
frequency techniques such as the Geometrical Theory of Diffraction 
(GTD), the Equivalent Current Method (ECM), the Aperture Integration 
(AI) and a modification of the Physical Theory of Diffraction ( PTD) . 

The edge diffraction, due to the rim of the open-end, is calculated via 
the ECM in which the equivalent currents are found from GTD. The AI 
together with a modified PTD approach is applied to get the modal 
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coupling into the interior and its subsequent radiation back into the 
exterior region. Since there are many propagating modes inside the 
cavity, it is expected that only a few modes which are closely phase 
matched to the direction of the incident wave are strongly excited. 
This claim is tested and comparisons with measurements are also shown. 

The experimental measurements are performed on a model at the 
compact range facility of the ElectroScience Laboratory of the Ohio 
State University. 
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CHAPTER II 


ANALYSIS OF ELECTROMAGNETIC BACKSCATTERING 
FROM AN INLET STRUCTURE 


SECTION I 



Figure 1. The geometry of the inlet model. 

The inlet model of interest in this work is shown in Figure 1. It 
is basically an open-ended cavity composed of two waveguide sections. 

The first section is part of a linearly tapered waveguide with one end 
open; the other end of this section is connected to the second section 
which is a uniform waveguide with a planar termination inside. The 
exterior of the second section is curved at the back end to minimize the 
scattering coming from the exterior features of the structure. The axis 
of the inlet coincides with the z-axis of the coordinate system, and the 
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y-axis is in the vertical direction. The model is made of wood and then 
made conductive using a conductive paint. The dimensions of of the 
inlet are shown in Figure 2. 

The side walls in each of the waveguide sections are parallel to 
the y-z plane; therefore, there is no tapering effect to be included for 
those walls. The first step in the analysis is finding the modal field 
expressions in rectangular and linearly tapered waveguide sections. 

This is done in Appendices I and II where the expressions are given, and 
the relationship between the mode sets of linearly tapered and 
rectangular waveguides is shown. Therefore, the effect of tapering 
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(b) TOP VIEW 


Figure 2. Side and top view of the inlet model. 


4 



can be included by modifying the propagation constant in the rectangular 
waveguide mode expressions which are transverse to the x-di recti on. So, 
the analysis is first done on a rectangular cavity as shown in Figure 3, 
and the solution is modified subsequently to include the effect of 
tapering. 

The geometry of the waveguide structure is shown in Figure 3. The 
cross-sectional dimensions of the waveguide are "a" and "b" in the x and 
y-coordinate directions respectively. The length of the waveguide from 
the open end to the back wall is given by the dimension L. All of the 
walls of the cavity are assumed to be perfectly conducting. The 
structure is illuminated by an incident plane wave (eJ wt time convention 
is assumed and suppressed) which is given by 


- < E > 




ikCxsine 1 cos (I) 1 


+ ysin^cosi )) 1 


+ zcose 1 ) 


( 1 ) 


where 0 < Q 1 < tt/ 2, 0 < <(>i < 2 tt are the elevation and aspect angles of 
the incident field direction, respectively, and k is the free space 
wavenumber. 


y 



Figure 3. Geometry of an open-ended rectangular waveguide cavity. 
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The scattered field is composed of two main contributions; one of 
these is due to the field scattered from the open end by itself, and the 
other is due to the field which is radiated from the open end. The 
latter undergoes multiple reflections between the open end and the 
termination at the back wall after it is initially coupled into the 
waveguide from the incident field. The fields resulting from these 
multiple interactions may be expressed as a convergent Neumann series as 
done by Pace and Mittra [1]. The same result is directly obtained by an 
alternative procedure based on a self-consistent method [2], However, 
from experimental measurements, these multiple interactions are 
determined to be negligibly small for the inlet model being considered 
here; therefore, their effect is ignored in this study. Only the first 
order interaction is discussed which contains the coupling of the 
incident energy into the interior waveguide modes through the open end, 
and subsequent reflection of these modes from the back wall and finally, 
the radiation of these reflected modes from the open end. In Appendix 
III, it is shown that the mechanisms of coupling into and radiation from 
the open-end are equivalent via the reciprocity principle. 

For a finite cross-section, the field incident at the open end 
excites a finite number of modes which propagate in the -z direction 
without attenuation, as well as an infinite number of evanescent modes 

A 

which attenuate exponentially away from the open end in the -z 
direction. In this study, it will be assumed that the dimensions "a" 
and "b" are large enough to excite at least a few propagating modes and 
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the length "L" is long enough for the effects of the evanescent modes to 
be negligible. 

In the third chapter of this report the analysis of direct edge 
diffraction from the rim of the open end is discussed; whereas, the 
discussion of the wave interaction between the open end and back wall 
is given in Chapter IV. 
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CHAPTER III 


DIRECT SCATTERING FROM THE RIM AT THE OPEN END 


In this section, the direct scattering from the rim (edge) at the 
open end is calculated using techniques based on Geometrical Theory of 
Diffraction (GTD). According to the GTD, the scattered field is 
initiated from some distinct points (diffraction points) on the rim edge 
as well as from the corners of the rim, respectively as a result of 
Keller's generalization of the Fermat's principle [3]. In addition to 
the singly edge and corner diffracted fields, there are multiply 
diffracted fields which are produced by rays that undergo multiple 
diffractions across the aperture. These multiple interactions may 
become important if the aperture dimensions are not sufficiently large 
in terms of the wavelength. 

The diffraction points migrate around the rim edge as the 
observation point changes position. In some cases, there may be a 
continuum of diffraction points contributing to the scattered field 
which produces a cuastic of the diffracted rays. In the case of 
backscattering when the incident field direction is parallel to x-z (or 
y-z) plane, every point contributes to edge diffraction. In order to 
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get a bounded result, one can use the so-called Equivalent Current 
Method (ECM) [4] based on the GTD. Away from the caustic regions where 
GTD is valid, the ECM generally blends into the GTD solution provided 
the inlet opening is sufficiently large in terms of the wavelength. The 
use of ECM also automatically, albeit approximately takes into account 
the presence of the corners at the waveguide opening. 

In the ECM, the equivalent currents I e q and M e q of the electric and 
magnetic type, respectively, are located at the rim, and they radiate in 
free space to give the diffracted field. The strengths of the 
equivalent currents are calculated from the GTD, but since they are 
incorporated in an integral, they give bounded results in the caustic 
regions of the GTD. 

The strengths of equivalent electric and magnetic currents are 
given by [4] 


I 


eq 


e • E 1 (edge) f^~ A 

Z Q sinB 1 . singj D s^’^ ) v jk e 


( 2 ) 


and 



e - H 1 (edge) /ST" - 

Y Q sing. sins^ D h^’^ ^ v j k e 


(3) 


where the half-plane diffraction coefficient is given by 


_ e -jir/4 

D s U,V) = 2 /2irk 
h 


cos 




hr 


cos 


(4) 
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which can also be written as follows: 

“sin ip/2 sin ip'/2“ 

_ e -jTf/4 cos ip/2 cos ip'/2_ 

^* rS ) = /2nk cos ip + cos ip' (5) 

Lh J 

The other parameters in Equations (2) and (3) are defined as follows: 

A 

e : unit vector along the edge direction 

Z 0 = Y 0 _1 : free-space impedance, and 

&i , (Bj) : the angles between the edge and incident 

(or diffracted) field direction. 

To define the angles ip 1 and ip, one first projects the vectors in the 
incident and diffracted field directions into a plane perpendicular to 
the edge direction. Then ip' is the angle between the half-plane surface 
and the projected incident field direction. The angle ip is similarly 
defined. The angles ip' and ip are measured from the inner faces of the 
waveguide. 

The radiated electric field is computed by 

jkZ Q e -jkR 

E Sl =TT / [R X R X 5 eq * Y 0 * x % 1 T - <f*' (6) 

rim 

— A 

where R = RR is the vector pointing toward the observation point from a 
source point on the rim. If the observation point is very far away from 
the open end, one can make the following approximations: 

A A A 

R s r - r • F' and R = r (7a;b) 
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where r=rr is the vector pointing toward the observation point from the 

A A 

origin. Also, r'=x'x+y'y is the vector from the origin to the source 
point. The unit vector in the observation direction is defined by 

^ A A A 

r = x sinecos<{> + ysinesin<j> + zcose (8) 

where 0 and <|> are the elevation and azimuth angles. 

As shown in Figure 2, the edges are numbered from 1 to 4. Singly 
diffracted fields from each edge are calculated separately, and total 
edge diffraction from the rim is found by adding the results for each 
edge. 



Figure 4. Numbering of the edges. 
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i .) Singly-diffracted field contribution from edges 1 and 2 

Edges 1 and 2 are the upper and lower horizontal edges as shown in 
Figure 4. Their contributions are similar in the far field; therefore, 
the expressions for edge 1 are given first, then the total contribution 
from both 1 and 2 are found. 

According to the above definitions, one finds that 

cose 1 cose 

cos ^‘ = - cos * = -i7^ ( 9a » b > 


sin ^i ^l-sin 2 eicos 2 <|>i sin0 d /i_sin 2 ecos 2 (|> (10a;b) 


7 (x',y')_ 
eq 


and 


[E’ e coseW+Ey-sW)] e Jk(x ’ sin6 ' cos » ,t > ,,sin0,sin *’ > 


Z„ sin 0. sine. 

Ola 


>»+') /jlT 


D $ ( ip, ) J ik x ( 11 ) 


[El coseW - e!(W)] eJMx'sIne’cos^ty 1 sine'sin*' ) 

<p U 


eq 


sing^ sin8 d 




IW~ 
.♦')V jk 


(12) 


Substituting Equations (11) and (12) into Equation (6) and adjusting the 
domain of integration for edge 1, one obtains that 
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,-jkr 


;si _ii. 


= 4irr 


sing^smgj 


8tt 

jk 


(EgCos 0 1 cos 4) 1 cos ecos 4>-si n 0 1 cos 0cos <|>E^) D g ( ij>, V ) 


+(-E^cos0 1 cos<j> 1 sin<}> + EgSin^sin^) D^(ip,t|>') 


+<f> 


(-Ejjxose’cos^sint)) - E^sin^ 1 si n 4 >) D_ ( v^, ip* ) 


+(-E^cos0 1 cos<|> 1 cosecos<|) + E qSI n <j> 1 cos 0cos <J>) ( ip, ip* ) 


i . .1 


j a e jk[x‘ (sinO^os^^singcosiji) + b( si no 1 sin <p^ +si n 0si n <p) ] dx , 


(13) 


This expression simplifies in the case of backscattering, namely 
0 1= 0 and <p n =<p - As mentioned above, the contribution from edge 2 is 
similar to Equation (13), the total contribution from edges 1 and 2 for 
the backscatter case is given by 
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For the backscattered field scan in the x-z plane, where d>i =0 , 
ip'=ir, one obtains that 

— ikr ^3 i p ( i 

_ . A e J j y ~ sine e sm(kasine ) 

^12 = " ® kr e tt -j • (15) 

sine 

Since both horizontal edges scatter always in phase in this plane, 
the singly diffracted field contribution of horizontal edges is 
independent of the dimension b. For an observation scan in the x-z 
plane none of the points along the horizontal edges satisfy the 
Generalized Fermat's principle. Therefore, there is no edge diffraction 
contribution from horizontal edges. But, since the edges are of finite 
length, the corners give a scattered field in every direction. At this 
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time there is an available approximate, uniform corner diffraction 
coefficient which is partly based on heuristic arguments [4]. If this 
corner diffraction coefficient is employed, one will simply obtain the 
result in Equation (15) as long as one is away from the shadow 
boundaries of edge diffracted fields. However, since the latter will be 
true in the present work, this ECM approach will yield the same result 
as that obtained from the use of the above corner diffraction 
coefficient. It is noted that in the far zone, only the corner 
diffraction effects dominate the total rim scattering; i.e, only the end 
points of the rim edges dominate rather than the whole length of the 
edges except in the caustic region where the entire length of the rim 
contributes to the field. 


i i . ) Singly-diffracted field contribution from edges 3 and 4 
In this case. Equations (2) and (3) take the following form: 


(eWsW + E'W) e jk(as1r,eW + y'sino'sW ) 

y cp 


eq 


Z Q sin6 i sine d 


/8tt * 

D-('lvP') v jk y 


(16) 


(E^cose’sin ^ 1 

<p 


E^costfi 1 ) eJ^asine^os^ 1 + y ' si n e 1 si n^ 1 ) 


eq 


sin^ sin6 d 


f 8tt 


D^( 'j'.'P' ) \J jk y 


(17) 
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where g. , 3 ^, and V are as defined in Equations (9) and (10). 

Substitution of Equations (16) and (17) into Equation (6) and 
defining the parameters for edge 3 one obtains 


£Sl 

b 3 


-jkr 

jk e 

^ irr sing.-singj 

i d 


/8tt 
V J‘ k 


(EgCose^in^cose sin<f> + E^cos^cosesin^) D s ( ip, ip * ) 


+ (E^cos^sin^coscj) - E gCOS <p 1 cos <p) ( ip, ip' ) 


i ... J 


+ <p 


(E^cose^in^^osef) + E^cos cp 1 cos <p) D ( \p, ip 1 ) 


+ (-Ejcose^in^cose si n <p + E^cose’cose si n <p) D h ( ip, ^ ' ) 


j b e jk[a(sine i cosc() i +sinecos<)>) + y^sine^in^+sinesin^)] ^ , 


(18) 
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The contribution of edge 4 is similar to Equation (18). The total 
contribution from edges 3 and 4 for the backscatter case is given by 


-si 


1 e "' , * <r 2cos(kasin9 1 cos<j>i ) sinCkbsine^sin^ ) 


E 34 “ sin 2 8. 


2cosij> 


si n 9 1 si n 4 ) 1 


ka 


kb 


e j 2 sinQ^os^ 1 j 2 si n Q 1 si n c}) 1 


■i 4 ijl 

E'cos 2 9 si n 2 4 > + E cos <t> si n 4 > cose )sin 2 2 


/r-i i i. i r-i ->i\ o ▼ 
•(E cose cos 4 > sin<|> - E p,cos z 4> ) cos z ^ 


+ <j> 


11111 1 ^ 
EgSi n 4 > cos9 + E^cos<j> ) cos sin 2 j 


1 111111 ^ 
+ (E^cos 2 0 si n<j> - E 0 cos4> sin<j> cos9 ) cos 2 j 


. (19) 


For an observation scan in the x-z plane, where <f> 1= 0, ^tt+9 1 , 8.- = 

I - 

one obtains 


gSl 

L 34 


e ‘ Jkr j sine 1 cosCkasIne 1 ) 


irr 


COS 9 


1 fl 1 

A A 0 ^ J ” 

<j>E^cos 2 2~ + 9 E 0 sin 2 Y~ 


( 20 ) 
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iii.) Double diffraction from vertical edges in the case of x-z 


plane scan. 


The singly diffracted rays from edges 3 or 4 will travel across the 
aperture and diffract once more at the other end of the aperture before 
going back to the receiver. Since in the x-z plane scan case, there is 
a continuum of doubly diffracted rays, the ECM is employed to obtain 
their contributions. 

Let E^, be the electric and magnetic fields in y-di recti on 
incident on the fourth edge due to the diffraction from third edge, then 
one obtains 

g-jka+jklsine 1 

E43 = y E* (at edge 3) D $ ( =ir+9 i , (21) 


and 

e -jka+jklsine n 

H 43 = - i V O E 1 0 (at edge 3) D h (* , =ir+9 1 , 7= . (22) 

This field E^^ is incident on edge 4, diffracts from edge 4 and 
constitutes a part of the total doubly diffracted field. For the 
backscatter case, these doubly diffracted fields have a caustic; 
therefore, one needs to utilize the ECM to obtain a bounded result. 
The equivalent electric and magnetic currents are given by 


-43 


D s (V=2> ^-Q 1 ) 


I 34 

eq 


/ 8tt 

V jk 


and 


(23) 


18 




(24) 



y • ^43 ° h (i|>'=2, ' }l=TT-0i ^ 

/ 8tt ~ 

y vlk y 


the doubly diffracted fields from these equivalent sources are given by 
the integral of Equation (6) such that 



The doubly diffracted fields due to diffractions at edges 4 and 3, 
respectively, are given by similar expressions. The total doubly 
diffracted field is found by adding these two contributions such that 


gS2 

b 34 


e -jkr be -jka 


irjk 2 cos 2 © 1 


, B 1 0 1 - 

y E^cos 2 2 ~ - 0E 0 sin 2 2 ~ 


. (26) 


In addition to doubly diffracted rays, there are multiply diffracted 
rays at the aperture. The analysis of multiple diffractions can be done 
in general by the GTD, provided that all multiply diffracted fields are 
ray optical. As is well known, the GTD diffracted fields are not ray 
optical in the transition regions [4]. In the case considered here one 
vertical edge is exactly on the reflection shadow boundary of the other 
edge, where it has to match the discontinuity of the GO (Geometrical 
Optics) field. Since it is exactly on the shadow boundary, one can 


19 



decompose the diffracted field into two ray-optical fields as done 
previously [5]. The subsequent (double) diffration of these two ray 
optical components of the singly diffracted fields can be done via 
conventional application of the GTD. The evaluation of the third order 
diffraction follows similarly and is given by; 



e-J kr be' Jk2a 
r 2ir /ujka 



cos 2_ 2 

cos 2 e i 


cosfkasinb 1 ) e j kasinel 



sin 2 2 ~ 
cos 2 © 1 


(27) 


The higher-order diffractions, beyond the third-order term, are not 
considered here because they are very weak in comparison to the lower 
order terms just analyzed. 


i v . ) Numerical results for Chapter III 

To check the validity of the analysis, the calculated results are 
compared with a set of measurements obtained using the Compact Range 
Facility at ElectroScience Laboratory, the Ohio State University. The 
measurements are done on an inlet model as shown in Figures 1 and 2. In 
order to remove the interior cavity effects, the inner surface of the 
back wall of the inlet is covered with absorbing material. Two sets of 
measurement results are obtained. 

a.) At a fixed frequency, the aspect angle ( 0 ) from the z-axis is 
varied in the horizontal (or vertical) plane from 0 to 90 degrees. 
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b.) At a fixed aspect angle in the horizontal (or vertical) plane 
the frequency response is measured, starting from 8 GHz up to 12 GHz at 
10 MHz steps. Therefore a bandlimited frequency response is obtained. 
Then, this response is processed through a Kaiser-Bessel window, and its 
inverse Fourier Transform is obtained using an FFT algorithm which 
represents a time domain response. Since the measured spectrum is 
bandlimited, the time domain response represents the impulse response of 
the target which is convolved with 
si ncu^t 

V t cosa) 0 t (28) 

where u> b =half -bandwidth (2 GHz) and u) q is the center frequency which is 
10 GHz in this case. 

The measured and calculated radar cross sections (RCS) versus the 
aspect angle are shown in Figures 5 through 16 for different 
polarizations and frequencies. In each figure measurements (solid line) 
and calculations (dashed line) are drawn on the same scale of dB over 
m 2 . Figures 5 through 10 correspond to a horizontal (x-z) scan; 
whereas. Figures 11 through 16 correspond to a vertical (y-z) scan. In 
both the horizontal and vertical scans, the radar cross section is 

A A 

measured and analyzed for two polarizations (<j> and e). Finally, each 
case is repeated at three different frequencies, namely around 8, 10 and 

A 

12 GHz. In both the horizontal and vertical scans, a <fi directed 
incident electric field will not diffract from two edges which are 
parallel to the scan plane. Therefore the contribution to the radar 
cross section in these cases is dominated by the diffracted fields from 
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RADAR CROSS SECTION (OB) 

-40. -30. -20. -10. o. 10 



THETfi (DEG) 


Figure 5. Aspect angle scan in horizontal (x-z) plane. Vertical U) 
polarization, frequency = 8.02 GHz. 

measured, — - calculated results. 






RADAR CROSS SECTION (DB) 



Figure 7. Aspect angle scan in horizontal (x-z) plane. Vertical (<j>) 
polarization, frequency = 11.95 GHz. 

measured, — calculated results. 
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THETfl (DEG) 


Figure 9. Aspect angle scan in horizontal (x-z) plane. Horizontal (e) 
polarization, frequency = 9.98 GHz. 

measured, — calculated results. 










THETfi (DEG) 


Figure 11. Aspect angle scan in vertical (y-z) plane. Horizontal (J) 
polarization, frequency = 8.02 GHz. 

measured, — calculated results. 







30. 60. 

THETA (DEG) 

Figure 13. Aspect angle scan in vertical (y-z) plane. Horizontal (J) 
polarization, frequency =11.95 GHz. 

measured, — calculated results. 





RADAR CROSS SECTION (DB) 


Ml/il 


30. 60. 

THETA (DEG) 

Figure 15. Aspect angle scan in vertical (y-z) plane, 
polarization, frequency = 9.98 GHz. 

measured, — calculated results. 



30. 60. 

THETfl (DEG) 

Figure 16. Aspect angle scan in vertical (y-z) plane. Vertical (8) 
polarization, frequency = 11.95 GHz. 

measured, — calculated results. 






the remaining two edges. As shown in Figures 5 through 7 and 11 through 
13, the calculations agree reasonably well with the measurements. 
However, for the other polarization (e), all four edges of the open end 
contribute to the radar cross section; however the contribution from two 
of these four edges reduces essentially to that from the end points (or 
corners). As shown in Figures 8 through 10 and 14 through 16, the 
agreement between the measured and calculated results is not as good as 

A 

the ^-polarization case, especially the horizontal scan. The reason for 
this discrepancy will become clearer after discussing the frequency 
domain responses, suffice it to say for now that it is due to the 
imperfections of the model. 

Finally one notes that the rim scattering analysis performed here 
is valid for aspect angles away from the shadow boundaries of the 
diffracted field. For aspect angles close to the shadow boundaries of 
the edges (e close to 90 degrees) the approach described here should be 
modified; this modification will not be given here. 

Frequency scans are done in the horizontal (<|>=0) plane, with two 
different polarizations and with the angle 9 fixed at four different 
values, namely 0, 15, 30 and 45 degrees. The measured (solid line) and 
calculated (dashed line) radar cross section results are plotted on the 
same graph and given in dB over m 2 . The phase variation of the radar 
cross section is also shown in the figures. To be able to make a phase 
comparison, the measured data has been processed so that the two results 
have the same phase center. The frequency spectra are then inverse 
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Fourier transformed to obtain the time domain responses. As explained 
before the time domain response is not an impulse response although it 
has been denoted like that on the plot for ease in understanding. 
Therefore, to make the comparison easier, the curves corresponding to 
measured and calculated time domain returns are shifted by equal amount 
from the center line. The scale in the time domain is dimensionless and 
should be taken as a relative scale. Finally, due to the malfunctions 
of electronic instrument during measurements, there are some glitches in 
the measured data. They are kept as they appeared originally, since 
their presence does not effect the overall characteristic of the curves. 

As in the case of aspect angle scans, the amplitude and phase of 
the measured and calculated RCS results are in good agreement especially 

A 

for a <j> polarized field (Figures 17 through 20). Also, it can be seen 
from the time domain plots that almost all of the measured return is 
from the rim at the open end. The absorber put at the back of the line 
very effectively removed any internal cavity effects. 

In the case of the horizontal (0) polarization, (Figures 21 through 
24) the time domain plots show a relatively large return occuring later 
in time than the return from the open end. The variation of the 
position of that return with different 0 values implies that it is 
coming from the external surfaces of the inlet model. Therefore as 
shown in Figures 22 through 24, the comparison between calculated and 
measured radar cross section results is not good. This also explains 
the discrepancy in the aspect angle scan results of Figures 8 through 10. 
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RADAR CROSS SECTION (DB) 

-WO. -30. -20. -10. 0. 1 



FREQUENCT IN GHZ 



Figure 17a. Frequency scan in horizontal (x-z) plane 

A 

<j>=0°, 9=0°, vertical (<j>) polarization. 
measured, — calculated results. 








Figure 18a. Frequency scan in horizontal (x-z) plane. 

4>=0° , 9=15°, vertical (<{>) polarization. 
measured, — calculated results. 
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IMPULSE RESPONSE 



+6. -14. -2. 0, 2. H. 6. - 8. 10. 12. 14. 


TIME IN NHNOSECS 


Figure 18b. Inverse Fourier transforms. 
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FREQUENCY IN GHZ 



Figure 19a. Frequency scan in horizontal (x-z) plane. 

<|>=0 o , 9=30°, vertical (<j>) polarization. 
measured, — calculated results. 
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Figure 19b. Inverse Fourier transforms. 


41 




' 8.0 9.0 10.0 11.0 12.0 


FREQUENCY IN GHZ 



FREQUENCY IN GHZ 


Figure 20a. Frequency scan in horizontal (x-z) plane. 

<j>=0°, 0=45°, vertical ($) polarization. 
measured, — calculated results. 
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IMPULSE RESPONSE 



Figure 20b. Inverse Fourier transforms. 
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Figure 21a. Frequency scan in horizontal (x-z) plane. 

<p=°°, 9=0° , horizontal ( 9 ) polarization. 
measured, - — calculated results. 
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Figure 21b. Inverse Fourier 
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Figure 22a. Frequency scan in horizontal (x-z) plane. 

4>=0°, e=15°, horizontal (e) polarization. 

measured, calculated, 

measured (time gated) results. 
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Figure 22b. Inverse Fourier transforms. 
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Figure 23b. Inverse Fourier transforms. 
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Figure 24a. Frequency scan in horizontal (x-z) plane. 

<t>=0° , 9=45°, horizontal (e) polarization. 

measured, calculated, 

measured (time gated) results. 
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To compare the returns coming from the open end only, the unwanted 
return is gated out from the measured time domain data and is indicated 
as such in each figure. The actual measured and calculated returns are 
shown with constant shifts, as before. These time domain plots were 
then converted into the frequency domain and shown on the same graph. 

As can be seen in Figures 22a, 23a, and 24a the calculated (long-dashed 
line) and gated-out measured (short-dashed line) results agree fairly 
well since they both correspond to the returns pertaining to the rim of 
the inlet. 

As explained before, the time domain responses shown are not the 
actual impulse responses because of the bandlimited nature of the data 
in frequency domain. If one had more frequency domain information then 
the time domain results would be closer to true impulse responses. To 
show this, calculations corresponding to Figure 20a are repeated and 
expressions of Equations (20), (26) and (27) are extended down to 10 
MHz. The resulting frequency domain plot is shown in Figure 25a. Then, 
this data is inverse Fourier transormed, and the time domain result is 
shown in Figure 25b. This result is closer to an actual impulse 
response, and it clearly shows the single, double and triple 
diffractions from the vertical edges of the open end. Note that the 
double order diffractions originating from both edges return to the 
receiver at the same time. 
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Figure 25b. Inverse Fourier transforms. 



CHAPTER IV 


THE ANALYSIS OF INTERIOR CAVITY EFFECTS 

As described at the beginning of this report, the coupling of 
incident energy into the interior waveguide modes and their subsequent 
radiation after undergoing multiple relfections between the back wall 
termination and the open end comprise the cavity effects. The 
relationship between the analysis of the problem of modal coupling into 
and radiation from the open end is given in Appendix III. The classical 
approach to determine the radiated field from the open end due to a 
given mode involves an integration of the modal field across the 
aperture. The edge effects inherent in this procedure are not as good 
as predicted by the GTD. However, the edge effects given by the 
aperture integration method can be improved by the modified PTD approach 
as described in Appendix IV. The underlying argument in Appendix IV is 
that one can improve the aperture integration by adding the Ufimtsev's 
equivalent currents at the edges. 

The modal field expressions in the rectangular waveguide as shown 
in Figure 1 is given in Appendix I (Equations A1-A13). Each mode is 
decomposed into four plane waves, each of which is a ray optical field. 

A 

The transverse (to z) components of the fields are as follows: 
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Nnm 

E = U 

modal 


a a 


-jn a x-jm b y-jS nn ,z 


(-xU+yV) e 

jn a x-jm b y-j6 nm z 


+ (-xU-yV) e 

. - J n a x+ jV j0 nm z 


+ (xU-yV) e 


+ (xU+yV) e 


-jn a x + jm b y-je nm z 


(29) 


and 


Ynm Nnm 

H = “43 — 

modal 


+ (xV-yU) e 


(-xV-yU) e 
j n a x-jm b y-jB nm z 
j" a x+jra b y-j B nm z 


-jn a x-jm b y-jB niI1 z 


+ (xV+yU) e 

. . -jnx+jm h y-j B z 

+ (-xV+yU) e a b nm 


(30) 


The directions of propagation of these four rays are given by the 
exponents. 

The integration of the modal fields across the aperture provides an 
approximate radiated field which is given by 


JkZ o a b 
E a,nn, 4 " J l 


RxR x j ”° dal + Y R x M ™ oda1 
eq o eq 


; -jkR 

“ R dy'dx' 

(31) 


56 


where 


5 model 1 = * h -N Y 

eq modal = mn nm 


xUcosn a xsi nm b y-yVsi nn a xcosm b y 


(32) 


and 


|jj modal 
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^modal xz 


-N, 


mn 


xVsinn a xcosmj 3 y+yllcosn a xsinm b y 


(33) 
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o eq 
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(34) 


In the far field 


R = r - r 1 • r 


and 


RxRxJ_ l ! K3da1+ YjxM™ dal = -N nm Ucosn a xsinm b y 


(35) 
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(36) 
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Substituting into Equation (31), one obtains that 

~r kN nm -ikr i ksine , 

E a,nm = " “Swr e e ~ 2 (acos<j>+bsin<fr) 


* , Y nm x » Y nm/ 

■ 0cos ^>( l+y^”cos 6 ]+())Si n i(i(cos 0 +- y^-j 
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j ksin 0 cos<j)+n a + ( “ J ) ksin 0 cos 4 >-n a 


' m si n 7 (ksin 0 co S(j)+mb) m sin 2 tksin 0 cos<|>-mb)~ 
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(37) 

. sinx 

As it is seen from Equation (37), the pattern has a type behaviour 

* * _ sinx 

in both 0 and <j> directions. There are four difference ~ * forms, each 

of which is due to the integration of a plane wave component of the mode 

, % sinx 

as ( described in Equation (29). Therefore each ~ — peaks up at the 
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corresponding plane wave directions. Also, as the waveguide dimensions 
get larger, the peaks become sharper. As a result, only a few modes, 
contribute strongly to the radiation around a given observation 
di recti on. 

As mentioned earlier, the edge effects given by aperture 
integration are improved using a modified PTD approach as discussed in 
Appendix IV. The equivalent electric and magnetic Ufimtsev edge 
currents are [6] given by 



A “ -j 

e • E (edge) 
Z Q sing.sing d 



and 



e • H 1 (edge) 
Y Q sing. sin& d 



(38) 


(39) 


where the Ufimtsev diffraction coefficients are 
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(40) 


after simplification, it can be written as follows: 
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(41) 
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Therefore, one obtains 
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(43) 


The radiated field for these edge correction terms is then given by 

jkZ o e _jkR 

^u,nm ■ TT / (R X R x i u t Y 0 R x M u ] -j- df . (44) 

rim 

In the far field, the above integrations for the radiated field can be 
done analytically in closed form and the result is; 
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( 45 ) 



The total radiated field is the summation of Equations (37) and (45); 
namely. 


= E r + E r 
a,nm u,nm 


(46) 


The total radiated field is the superposition of the radiation of all 
propagating modes. 


Numerical Results for Chapter 


The experimental measurements and calculations are done in two 
categories as in Chapter III; namely, aspect angle and frequeny scans 
and they are given in dB over m 2 . Since it is difficult to 
experimentally isolate the open end rim scattering from the cavity 
effects, the calculations include both effects for comparison purposes. 

The waveguide model is large in terms of the wavelength; therefore, 
a large number of propagating modes can exist inside the cavity. For 
example, at the frequency of 10 GHz, there are 172 propagating modes as 
well as an infinite number of evanescent modes. Since the waveguide 
axial length is very long with respect to the wavelength, the evanescent 
modes will decay very significantly such that they can be neglected. 

Then the next question is whether the propagating modes have a preferred 
direction of radiation in that one can choose a few special ones rather 
than all the propagating modes and include only those in the analysis. 

It is obvious from Equation (47) that modal radiation occur at discrete 
plane wave directions for each mode. Since each mode has a discrete 
radiation direction, some will, radiate close to desired direction, and 
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others will not. This being the case, one can anticipate that only a 
few inodes are significantly excited by the incident plane wave and 
radiate significantly in the backscatter direction. This claim has been 
checked numerically, and the results are shown in Figure 26. An 
open-ended rectangular waveguide with dimensions equal to the 
experimental model used in this study is analyzed at 10 GHz. The 
backscatter field is calculated in (x-z) plane by varying the aspect 
angle 0. Note that only the modal effects are included and the incident 

A 

field is assumed to be ^-polarized. It is found that only the TE „ 

no 

modes are excited in this plane for this polarization. Their plane wave 
directions or modal ray angles are tabulated in the figure. The 
scattered field is calculated by including all modal contributions as 
indicated by the solid line. For comparison purposes the contributions 
of TE_ (n=l,...5) modes and TE„„ (n=6...13) modes are shown in the 
same figure. Note that the first five modes radiate strongly and almost 
replicate the solid curve in the region close to their modal ray angle 
directions. The remaining seven modes are major contributors for the 
large aspect angles. After observing this fact, the research focused on 
finding a rule of thumb procedure which could be used to select the 
minimum number of modes needed for a given direction of incident plane 
wave. From numerical results and comparisons with the experimental 
data, it was found that in horizontal and vertical planes only three 
modes were sufficient for the inlet under test. These three modes are 
selected based on their modal radiation direction such that they are 
closest to the incident plane wave direction. The aspect angle scan 
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Figure 26. Importance of modes whose modal ray angles are near the 
angle of incidence. 
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results are shown in Figures 27 through 38. As done in Chapter III, the 
experimental measurement results are indicated by the solid line and 
calculations by the dashed one. In each case a comparison is provided 
for the calculations which include all modes versus three modes. In all 
cases a 1 dB/bounce energy loss is assumed to model the imperfection of 
the conductivity of the model used in the measurements. 

Since this new concept proved to be so valuable in the principal 
plane, it is next applied to <f>=45° plane to see if it fails when the 
incident plane wave direction is not aligned with the structural 
symmetry. The aspect angle scan is calculated in the <j>=45° plane for 

A A 

both the <j> and e polarized incident fields. The results are shown in 
Figures 39 and 40. In this case, the 18 preselected modes are compared 
with the complete 172 modes as shown in each figure. Note that more 
terms are necessary in this general case as indicated by the results 
shown in Figure 41 where the six mode result is compared with the 18 
mode one. Even so it is clear that one can use far fewer than the 
complete set of modes. 

Finally, the frequency scan results are calculated and compared 
with measurements in Figures 42 through 49. As in Chapter III, the time 
domain results are deliberately shifted to better illustrate the 
comparisons. Note that in each case the comparisons are exceptional. 

In Figure 50, the contribution of each mode is plotted separately 
for the geometry associated with the results of Figure 45. It is seen 
that as the mode number increases the given mode undergoes more bounces 
inside the waveguide, travels a longer distance and therefore its return 
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THETfl (DEG) 


Figure 27a. Aspect angle scan in horizontal (x-z) plane. 

Vertical (<j>) polarization, frequency 8.0 GHz 
measured, — calculated results. 








Figure 27b. 


Calculated aspect angle scans corresponding to Figure 27 . 

All modes are included 

xxx Only 3 modes are included 
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' THETfl (DEG) 

Figure 28a. Aspect angle scan in horizontal (x-z) 

a 

Vertical ($) polarization, frequency 
measured, — calculated results. 
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plane. 
10.0 GHz. 



THE T fl (DEG) 


Figure 28b. Calculated aspect angle scans corresponding to Figure 28a 

All modes are included 

xxx Only 3 modes are included 






THETfl (DEG) 


Figure 29a. Aspect angle scan in horizontal (x-z) plane. 

Vertical (<j>) polarization, frequency 12.0 GHz 
measured, — calculated results. 








Figure 30a. Aspect angle scan in horizontal (x-z) plane. 

* 

Horizontal (9) polarization, frequency 8.02 GHz. 
measured, — calculated results. 
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RADAR CROSS SECTION (DB) 



THETR (DEG) 

Figure 30b. Calculated aspect angle scans corresponding to Figure 30a. 

All modes are included 

xxx Only 3 modes are included 
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RflDfiR CROSS SECTION (DB) 



Figure 31a. Aspect angle scan in horizontal (x-z) plane. 

a 

Horizontal (8) polarization, frequency 9.98 GHz. 
measured, — calculated results. 
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THETfl (DEG) 


Figure 31b. Calculated aspect angle scans corresponding to Figure 31a. 

All modes are included 

xxx Only 3 modes are included 
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THETA (DEG) 

Figure 32a. Aspect angle scan in horizontal (x-z) plane. 

A 

Horizontal (e) polarization, frequency 11.95 GHz. 
measured, — calculated results. 
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o' ' ' ' 

inO. 30. 60. 90. 

* THETA (DEG) 

Figure 33a. Aspect angle scan in vertical (y-z) plane. 

A 

Vertical (9) polarization, frequency 8.02 GHz. 
measured, — calculated results. 
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Figure 33b. Calculated aspect angle scans corresponding to Figure 33a. 

All modes are included 

xxx Only 3 modes are included 
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RADAR CROSS SECTION (DB) 



Figure 34a. Aspect angle scan in vertical (y-z) plane. 

A 

Vertical (e) polarization, frequency 10.01 GHz. 
measured, — calculated results. 
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RfiDRR CROSS SECTION 
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Figure 34b. Calculated aspect angle scans corresponding to Figure 34a. 

All modes are included 

xxx Only 3 modes are included 
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Figure 35a. Aspect angle scan in vertical (y-z) plane. 

A 

Vertical (e) polarization, frequency 11.96 GHz. 
measured, — calculated results. 
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THETR (DEG) 


Figure 35b. Calculated aspect angle scans corresponding to Figure 35a 

All modes are included 

xxx Only 3 modes are included 





Figure 36a. Aspect angle scan in vertical (y-z) plane. 

A 

Horizontal ($) polarization, frequency 8.02 GHz. 
measured, — calculated results. 
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uO. 30. 60. 

' THETR (DEG) 

Figure 36b. Calculated aspect angle scans corresponding to Figure 36a. 

All modes are included 

xxx Only 3 modes are included 
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i^O. 30. 60. 

' THETfl (DEG) 

Figure 37a. Aspect angle scan in vertical (y-z) plane. 

Horizontal ( <f>) polarization, frequency 10.01 GHz. 
measured, — calculated results. 
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RRDAR CROSS SECTION (DB) 



Figure 37b. Calculated aspect angle scans corresponding to Figure 37a. 

All modes are included 

xxx Only 3 modes are included 
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90. 


is’O. 30. 60. 

THETA (DEG) 

Figure 38a. Aspect angle scan in vertical (y-z) plane. 

Horizontal (<(,} polarization, frequency 11.96 GHz. 
measured, — calculated results. 
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THETA (DEG) 


Figure 38b. Calculated aspect angle scans corresponding to Figure 38a 

All modes are included 

xxx Only 3 modes are included 




RADAR CROSS SECTION (DB) 



Figure 39. Calculated aspect angle scans in the plane <j>=45°, 

$ polarization, frequency = 10.0 GHz. 

All 172 modes are included (86 TE, 86 TM modes) 

xxx Only 18 modes are included (9 TE, 9 TM modes) 
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Figure 40. Calculated aspect angle scans in the plane <fr=45°, 

9 polarization, frequency = 10.0 GHz. 

All 172 inodes are included (86 TE, 86 TM modes) 

xxx Only 18 modes are included (9 TE, 9 TM modes) 
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Figure 41. Calculated aspect angle scans in the plane <j>=45°, 

9 polarization, frequency = 10 GHz. 

All 172 modes are included (86 TE, 86 TM modes) 

xxx Only 6 modes are included (3 TE, 3 TM modes) 
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Figure 42b. Inverse Fourier transform. 
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Figure 43a. Frequency scan in horizontal (x-z) plane. 

A 

<j>=0°, 9=15°, vertical (<{>) polarization. 
measured, — calculated results. 
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IMPULSE RESPONSE 



Figure 43b. Inverse Fourier transforms. 
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Figure 44a. Frequency scan in horizontal (x-z) plane. 

4=0°, 9=30°, vertical ({) polarization. 
— measured, — calculated results. 
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Figure 44b. Inverse Fourier transforms. 
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me.O 9.0 10.0 11.0 12.0 

FREQUENCY IN GHZ 


Figure 45a. Frequency scan in horizontal (x-z) plane. 

<j>=0°, 9=45°, vertical (<j>) polarization. 
measured, — calculated results. 
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Figure 45b. Inverse Fourier transforms 
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Figure 46a. Frequency scan in horizontal (x-z) plane. 

A 

<j>=0°, 0=0°, horizontal (e) polarization. 
measured, — calculated results. 
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Figure 46b. Inverse Fourier transforms. 
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Figure 47a. Frequency scan in horizontal (x-z) plane. 

a 

<J>=0°, 0=15°, horizontal (9) polarization. 
measured, — calculated results. 
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Figure 47b. Inverse Fourier transforms. 
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Figure 48a. Frequency scan in horizontal (x-z) plane. 

$=0°, 9=30°, horizontal ( e) polarization. 
measured, — calculated results. 


105 


IMPULSE RESPONSE 







juil 


Z 

m 

ASUR 

ED 






iff 












o 

> 

r* 

CULATED 






■ 




: 


w6. -4. - 2 . 0. 2 . 4. 6. 8. 10. 12 . 14. 

TIME IN NRNOSECS 


Figure 48b. Inverse Fourier transforms. 
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Figure 49a. Frequency scan in horizontal (x-z) plane. 

<j>=0 o , 9=45° , horizontal (e) polarization. 
measured, — calculated results. 
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Figure 49b. Inverse Fourier transforms. 
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IMPULSE RESPONSE 



Figure 50. The separate returns of modes for the case in Figure 45b. 
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is received later in time, causing a stronger modal dispersion. Also 
note that the modal terms add up to form ripples in the total result 

from the cavity effect. If not understood correctly one might think 

» • • • , 

this is indicating different scattering centers, which would not be 
true. This shows that modal propagation should carefully be traced if 
one wishes to obtain the true response for the cavity effects. 
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CHAPTER V 


CONCLUSIONS 

The goal of this research is to understand the dominant scattering 
mechanisms associated with a jet inlet structure and to develop some 
simple and accurate techniques to analyze them. Basically the 
backscattered field is composed of two major mechanisms; namely, the 
diffraction from the open end and radiation from the open end due to 
interior reflections or cavity effects. Both mechanisms are separately 
discussed, analyzed, and compared with aspect angle and frequency scan 
measured results. The analysis and numerical results for the 
diffraction by just the open end is presented in Chapter III. The 
results indicate that the GTD provides an accurate estimation of this 
diffraction term. Chapter IV deals with the analysis of interior cavity 
effects. The latter analysis is based on the expansion of the field 
inside the inlet in terms of rectangular waveguide modes. The coupling 
of the incident field into the waveguide modes and their subsequent 
radiation from the open end are obtained by an integration over the 
aperture at the open end. In appendix IV, it is shown that this 
aperture integration can be corrected in the same way as the P 0 is 
corrected via the PTD. However, the integration process does not give 
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any physical insight to the problem, and as the number of propagating 
modes increase with frequency it becomes cumbersome and inefficient. 
Therefore, in this research, new approaches to substantially improve the 
efficienty of these aperture integrals was investigated. It was shown 
that every mode can be decomposed into four ray-optical fields each of 
which is propagating in a particular direction. These directions are 
called the modal ray directions. It was determined that the modal 
radiation is strongest in these modal ray directions. Therefore, for a 
given radiation direction one can include only the radiation of those 
modes whose modal ray directions are closest to the radiation direction. 
Numerical results are included to indicate that the contribution of only 
three modes with modal ray directions closest to the given observation 
direction accurately approximates the contribution of all the 
propagating modes. This is an important result, because it combines 
modal radiation with ray-optics, and therefore it can be applied to many 
different and complex waveguide geometries to effectively select only 
the few significant modes from the entire set of propagating modes. 

Also, the same modal ray approach is useful in calculating the modal 
reflection and transmission from interior inlet discontinuities. 

The model studied here had a small interior discontinuity, and, as 
expected, it did not seem to significantly influence the final results. 
It is proposed as a part of future research, to build and analyze a 
model with more pronounced interior discontinuities which might occur in 
practice. 
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APPENDIX I 


MODAL FIELD EXPRESSIONS IN A RECTANGULAR WAVEGUIDE 



Figure Al. Rectangular waveguide geometry. 


As described in [7], the modes in a rectangular waveguide as shown 
in Figure A.l, can be classified into sets of transverse fields to a 
coordinate direction. 

A 

The mode sets transverse to the axial (z) direction are important 
and widely used, because it applies to uniform nonrectangular 
cross-section guides. However, in many problems, mode sets transverse 

A A 

to x or y coordinate may be more suitable. 
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These sets can be calculated from an axial or transverse vector 

A 

potential [7]. The TE nm and TM pm to z type modal field expressions are 
given by: 


TE Modes 


TE Modes 


N nm 3F ("H' cosn a xcosm b y 
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h 6 

where the normalization factors N and are defined so as to 

nm nm 


a b 

/ / e x R • z dx dy = 1 

o o 


(A9) 


and 


on 


1 if n = 0 

2 if n * 0 


(A10) 


n = 


nu 

» m b = b 


(All) 


Also, S nm = is the propagation constant in 

z-di recti on, namely; 


A A 


E modal = t e x x + e y y + e z z ] e 


+j 3nm z 


(A12) 


A A 


"modal - t\ x + h y y + h z Z 1 e 


+j enm 


(A13) 


Finally, Y 0 =Z~* is the free-space admittance. 
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The modal sets transverse 

to x-axis can be 

written as follows: 



T^x,nm Modes 

TMx.nm 

Modes 
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n a cosn a xcosm^y 


From Equations (Al) through (A19), one can easily conclude that 
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or, alternatively 
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APPENDIX II 


MODAL FIELD EXPRESSIONS IN A LINEARLY TAPERED WAVEGUIDE 

y 



Figure A2. Geometry of a linearly tapered waveguide. 


In this problem, it is possible to find mode sets transverse to 

A 

x-axis. The fields can be computed from the scalar function 


a 2 3 2 
3y 2 + 3z 2 + 



ip=0 


(A24) 


with appropriate boundary conditions, t is the x-component of the 
magnetic vector potential A in the TM X case and the x-component of the 
electric vector potential F in the TE X case [7], 
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Case 1 TE modes: These modes satisfy the relation E =0 and can be 

. —i— i. i. _ A A 


derived from F=xi|>. 




The field components are given 

by; 



-1 jhfc 
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The appropriate scalar function \J> is given by; 
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where 


2 .OTT 

k - (t) 


k.. = k - f— 


(A29) 


The explicit expressions for the outgoing field components are given by; 
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1 , mr ni imr d (?) t 

H p = jkZ 0 (a ) cos (a - *) cos C ) dp ( k tP) 
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(A34) 


(A35) 


For small (f> 0 and large p, the Debye Asymptotic approximation for the 
Hankel function [8] 
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can be used. For small z ; cos- 1 z -*■ ?, so the above relation becomes; 
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which goes to the large argument form [7] for z>>v. 
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If one makes use of Equations (A37) and (A38) in the field 
expressions of Equations (A30) through (A35), one obtains the following 
expressions (up to a common factor) for an outgoing wave; 
mu - 1 


E * P*o 'W s1n (5TX) 
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B nm p 
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(A41) 

(A42) 

(A43) 

(A44) 

(A45) 


is the function determining the phase advance in p-di recti on. Note that 
B nm is a function of p. By comparing Equations (A39) through (A44) with 
Equations (A14) through (A19), one can see that the expressions of 
slightly tapered waveguide fields may be approximated by TE X , TM x -type 
rectangular waveguide mode sets. To do this, one has to make the 
following approximations. 
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+ z 


(A46) 


P< 1 > + y (A47) 

p4> q -*■ b (A48) 

together with Equation (A45) where p-dependence cannot be approximated 
by Equation (A48), since it is the function determining the phase 
advance which is more sensitive to the approximations. 


Case 2 TM modes: These modes satisfy the relation H =0 and can be 

~ ~ — A ^ 

— A 

derived from A=xi|>. 


The field components are given by: 
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The appropriate scalar function is given by: 
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The explicit expressions for the outgoing field components are given 


by: 
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Again using Debye's Asymptotic form for the Hankel function, one obtains 
the following field expressions: 
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These modal expressions under the approximations defined by 
Equations (A46) through (A48) become equivalent to the TM X modes in a 
rectangular waveguide. 
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APPENDIX III 


APPLICATION OF THE RECIPROCITY THEOREM TO FIND THE COUPLING 
OF INCIDENT ENERGY INTO THE MODAL FIELD BY KNOWING THE MODAL 
RADIATION FROM THE OPEN END OF A RECTANGULAR WAVEGUIDE 



Figure A3. The geometry of the problem. 


Figure A3 shows an open ended rectangular waveguide geometry with 
perfectly conducting walls. Let there be modal fields inside the 

a 

waveguide propagating in +z direction and represented by: 


E* = c* (e +e ) e n 
n n ' n zn' 

-+ + ,, - , - j8 n z 

H n “ c n < fi n +(! Z n> e 


(A64) 

(A65) 
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where e n (or h n ) is the transverse electric (or magnetic) field, and i zn 
(or fi ) is the longitudinal electric (or magnetic) field, c* is the 
coefficient of the nth mode. 

Some of the energy carried by this modal field will be radiated into 
the free space through the aperture. The far field radiation can be 
calculated approximately by the aperture integration of the modal field 
and the resulting radiated fields are shown by E r and H r . 

The rest of the energy of the incident mode will reflect back to the 
guide in terms of an mth mode. 




(A66) 


. J ’ 0 m z 

H ‘ * r c! (-fL+fLJ e m 

m mn n m zm 


(A67) 


where r nm is the complete modal reflection coefficient from the open 
end. 

First consider the geometric plane defined by z=-L inside the 
semi -infinite pipe, and let S Q denote the area of this waveguide cross 
section at z=-L as shown in Figure A3. Then let S denote the surface 
area which tightly encapsulates the complete outer (exterior) surface of 
the semi -infinite waveguide, and also a portion of the inner walls of 
this waveguide up to the distance z=-L within the guide. The surface S 
does not include the plane at z=-L within the pipe. Let z denote the 
sphere at infinity which surrounds the semi -infinite guide such that S 
on the exterior or outer wall of the semi -infinite guide is connected to 
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the surface E at z+ -» as in Figure A3. Next, consider the following 

two cases. In the first case, the semi -infinite rectangular waveguide 

is excited from within by the modal fields with amplitude c* which 

* 

propagate in the +z direction. However, in the second case, the 
geometry is excited by external fields, E e and H e , which for convenience 
is assumed to be produced by an electric test source J g at Rp exterior 
to and in the far zone of the semi -infinite guide. 

Let the equivalent sources and M g be located at z=-L in the 
guide. 


J s = n 




.j»n L 


+ e r. 


m 


mn n 




-J'V-' 


(A68) 


= -n x 


j0 n L + _ _ 

c (e +e_ n ) e + E r c (e-e) e 
n ' n zn' _ mn n ' m zm' 

m 


( A69 ) 


where n is the unit vector pointing into the volume V, enclosed by 
surfaces S+S q +E. Form reciprocity one obtains: 


I I 

E+S+S, 


E r x H e - E e x H r 


nds = / / / J • E r dv 
V e 


(A70) 


However, by the boundary conditions 


n x E' i=0 


on S 


and 


(A71) 


n x E j=0 
on S 


(A72) 
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Also, (E r , H r ) and (E e , H e ) satisfy the radiation condition on E. 
Therefore Equation (A70) reduces to: 

- 11 H 8 

S 

o 

since 

= P e 6 * (*™) 

a point source in the far field. 

In Equation (A74) R p is the position vector from the reference point 
to point P. 

Inside the waveguide the fields E e and H e will have the following 
representation 

E ' e ■ l A P < e 'p - 5 zp> eJ6pZ < A75 ) 

H 6 = E A" (-R p + R zp ) e J6pZ (A76) 

Substituting Equations (A75) and (A76) into Equations (A73) one obtains: 


• M $ ds + /I E e . J s ds = p e • E r (R p ) , (A73) 

S o 


-J B n L . 

11 E A (-fi +fi ) e p c 

( n P P zp n 

5 o p 




e xz e + E r „(e m xz)e 
n „ mn v m ' 

m 




ds 


- IJ l A p <VV e 


■J®d L + 


J 3 n L » 

(z x R )e - E r (zxR )e 
n' mn nr 

m 


ds 


" p e • ? (V 


( A77 ) 
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Using the orthogonality properties of waveguide modes one 


|| e„ x fi • z ds 
n n 

b o 


(-c+ A") 2 = p • E r (IT) 
n n' r e p 


-2 c+ A" 
n ,,n 


P e • E r (R p ) 


// e p x fi n • z ds 

0 


Ijl order to produce an incident field E 1 
the waveguide, one can let 

i ^ i 4ir 

p e - - <0 E e * * E i> 1ST % e P 


A i A i 

9 E e + 9 E * at 


Thus, 


A -• A J ~p 

(6 E ; + ♦ Ej) • E r (R n ) R„ e p 


P' P 


A_ = 


+ — A 
-2 c || e n x fi • z ds 
n J n n 

*0 


4tt 

juy 


obtains: 

(A78) 


(A79) 


the opening of 
(A80) 


(A81) 
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APPENDIX IV 


CORRECTION FOR THE APERTURE INTEGRATION 



Figure A4. Scattering from a perfectly-conducting half-plane. 


As shown in the figure, a perfectly conducting half-plane is located 

in the region y=0, x<0. The half-plane is illuminated by a plane wave 
u 1 =e jkecos(*-*') (n/2 < 

Let S t , denote the top and bottom surfaces of the half plane, 
respectively and denote the aperture surface at x=0, y>0. 

a 

Since the geometry is infinite in z-di recti on the problem can be 
reduced to a scalar problem where the total field u may represent either 

A A 

the z-component of the electric field (E z ) or the z-component of the 
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magnetic field (H z ). The former case is called the soft case and the 
latter is called the hard case. The analysis for both cases are very 
similar, therefore only the soft case will be considered here. In this 
case, the total field E z can be written as the combination of three 
contributions; namely, the incident, the reflected and the diffracted 
fields. The incident and reflected fields are also called as the 
geometrical optics (GO) fields. The GO field is given by the following 
expressions: 


go _ _ j k pcos U+$ ') _ g j k pcos ($+*') _ e ^ k x x ?Jsin|( 1 y . Qn ^ ^ 


E- e 


= 0 


; on S b 


since is in the shadow region of the half-plane. In Equation 


(A82) 
( A82 ) , 


k x = kcos^' 


and 


k^ = ksin<j>' 


(A83a;b) 


The transverse fields can be calculated using 

~jkZ 0 H = (-y aJ- + x gj) E z (A84) 


Therefore, the transverse geometrical optics fields can be obtained 
using Equations (A82) and (A83). 


qo i J' k x x i 

' jkZ o H x = jk y e 2cosk / 


(A85) 


' J ’ kZ o H v° = " Jk x ^ X 2sink y^ 


(A86) 
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In the physical optics (PO) approximation, one uses the following 
radiation integral to get the scattered field using the GO fields on the 
illuminated side of the half-plane (S t ). 


E =* / ds ' 
S 


N kz 0 G 0 (nxH) + (nxE) x V G Q ] 


(A87) 


where 

G 0 =iH< 2) ( k | p-d ‘ | ) 


(A88) 


is the free-space Green's function in two-dimensions. 

A — 

One notes that on S { ; n x E = 0. Therefore the PO approximation to 
scattered field is given by: 

E^° = / ds' [-jkZ Q G q (yxH 90 )] ( A89 ) 

S t 

Using Equations (A85) and (A86) in Equations (A89), one obtains: 

00 

E z = -I [- jkZ 0 H x • G ol dX ' (A90) 


The Green's function Go can also be written as: 


-J k »y-j k v x 


G o ■ 


- TF J «k y -T7 


jKy‘ 


-jk x x-jkv 

j - « y J k ,»' 

e y = - / dk x ] j- e 

—oo * 

(A91) 


where 


k 


2 

x 



(A92) 
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Substituting Equations (A85) and (A91) into Equations (A90) 
changing the order of integration, one obtains: 

-Jk x-jk v 


-P0 zi 
'z 


x ° y- 1 


E-”=2¥ / dk 


i 1 r x 
k v / e 

J o 


jk^x* + jkx 1 


dx* 


-P0 lA 


oo e 


-jk x x -o k y y k i 


2 ? I 


“j dk x 

k +k 
x x 


Using the transformation 


k =kcosa k =ksina 

x y 


Equation (A94) becomes 

F P0 _ f -jkpcos(a-<t>) s1n *‘ — r d 

b 7 - 2tt J e COSa+COSd) dot 


where the integration path c a is shown below: 



Rea 


Figure A5. The integration path c^. 


and 


(A93) 


(A94) 


(A95a;b) 


(A96) 
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On the other hand, if one integrates the GO fields over and to 
get the scattered field, one obtains: 


E* = J ds' [-jkZ G q (xxy Hy 0 ) + (xxz if) x 7G Q ] (A97) 

Sa 


since GO field is zero on S^. 

Using Equations (A82), (A86) and (A91) in Equations (A97), one 
obtains: 


J 

“ Tir 


/ dk 


-jk y y-jk x x 


k +k 

- y y 


k -k 

y y- 


(k x- k x } 


(A98) 


Using the transformations in Equation (A95), one finally obtains: 

E k _ _ J. f -jk P cos(a-<j)) si "»' (Agq) 

2tt J e cosa+cosd da 


Equations (A96) and (A99) are exactly the same, therefore the Ufimtsev 
edge correction [6] to Equation (A96) can be used without modification 
to correct for the aperture integration result of Equation (A99). 
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